Error estimates on Nuclear Binding Energies from Nucleon-Nucleon uncertainties 
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Despite great theoretical efforts the NN interaction can only be determined with a finite precision, implying 
an error upper bound for nuclear masses. We analyze for the first time the problem of estimating the systematic 
errors related to the form of the potential and their impact on nuclear binding. To this end we exploit the concept 
of coarse grained interactions to typical nuclear wavelengths. Our estimate gives an error AB /A ~ 0.1 — 0.4MeV 
for the binding energy per particle and paves the way for ab initio calculations tailored to such a precision. 
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Since the early days of Nuclear Physics the NN interaction 
has played a major role in the description of the properties 
of finite nuclei. While abundant sets of np and pp data have 
been collected along the years and accurate theoretical anal- 
yses have been carried out since the mid-nineties the 
impact of the NN uncertainties in the Nuclear Many Body 
Problem remains an open challenge. This is of utmost im- 
portance as it determines a priori a lower bound on the inac- 
curacy of first principles calculations [6] and might help to ad- 
vantageously optimize the computational cost. While Nuclear 
binding energies are experimentally known to high accuracy 
AB = 0.01 — lOKeV, liquid-drop model inspired mass fit for- 
mulae yield a lower theoretical accuracy AB = 0.6MeV (see 
e.g. Refs. @,[H and references therein). In the present work 
we face the problem squarely from the NN side by deducing 
and propagating two-body systematic errors to provide a first 
theoretical a priori estimate of binding energy uncertainties. 

Error analysis of NN phase-shifts for several partial waves 
became first possible when the Nijmegen group [1] carried 
out a Partial Wave Analysis (PWA) fitting about 4000 exper- 
imental np and pp data (after rejecting further 1000 of 3(7- 
mutually inconsistent data) with ^ 2 /dof ~ 1. The fit fixed 
the form of the potential to be an energy dependent square 
well located at a distance of 1.4fm, a One-Pion-Exchange 
(OPE) and Charge-Dependent (CD) contribution starting at 
1 .4fm and a One-Boson-Exchange (OBE) piece operating be- 
low 2 — 2.5fm. Unfortunately, the required energy dependence 
becomes messy for Nuclear Structure calculations. At present 
there are a variety of NN (energy independent) potentials fit- 
ting a large body of scattering data with ^ 2 /dof ~ 1 |Q]-|3|, 
but surprissingly error estimates on potential parameters are 
not given. Whereas all these modern potentials share the ven- 
erable and local OPE and CD tail and include electromag- 
netic effects, the unknown short range components of these 
potentials display a variety of forms and shapes, local poten- 
tials yj], or nonlocal ones implementing angular momentum 
dependence O], energy dependence [ 1 ] or linear momentum 
dependence yl |4|, |5|] . While in principle p—, L— and E— non- 
localities are on-shell equivalent (see e.g. Ref. BSD for a proof 
in a 1 /Mn expansion) they reflect truely different physical ef- 



fects and generally one should consider them as independent 
quantities; any specific choice is biased and hence becomes a 
source of systematic errors. 

We distinguish as usual in error analyses two sources of 
uncertainties: statistical errors stemming from the data un- 
certainties for a fixed form of the potential, and systematic 
errors arising from the different most-likely forms of the po- 
tentials. Clearly, the total uncertainty corresponds to adding 
both in quadrature. In what follows it is adventageous to take 
the viewpoint of considering any of the different potentials as 
an independent but possibly biased way to measure the scat- 
tering amplitudes and/or phase-shifts. Because the biases in- 
troduced in all single potential are independent on each other, 
a randomization of systematic errors makes sense. Thus, the 
overall spread between the various phenomenological models 
with ^ 2 /dof ~ 1 provides the scale of the uncertainty. 

In Fig. Q] we show the absolute (mean-square) errors for np 
partial wave phase shifts due to the different potentials fit- 
ting scattering data with ^ 2 /dof ~ 1 fl-Ht] as a function of 
the LAB energy. As one naturally expects the uncertainties 
grow with energy and decrease with the relative angular mo- 
mentum which semiclassically corresponds to probing an im- 
pact parameter b ~ (L+ l/2)/p, with p = the 
CM momentum, making peripheral waves to be mostly deter- 
mined from OPE. These analyses stop at the pion production 
threshold so that one probes distances till b m [ n ~ 1/A = 0.5fm 
with A = ^Jm n Ms. Generally, the PWA statistical errors [1] 
turn out to be smaller than the systematic bands displayed in 
Fig. Q] This counter-intuitive result relies not only on the spe- 
cific forms of potentials which treat the mid- and short-range 
behaviour of the interaction differently but also on the fact that 
the fits are mainly done to scattering amplitudes rather than to 
the phase-shifts themselves. Our purpose is to quantify the 
impact of uncertainties in Fig.Q]on Nuclear Binding energies. 

The most direct way of analyzing binding energy uncer- 
tainties from randomized systematic errors would be to un- 
dertake large scale ab initio calculations using the differ- 
ent forms of the set of N two-body potentials, say V 2 with 
i = 1, . . .N, yielding Z?W (A) whence a mean B(A) and a stan- 
dard deviation AB(A) can be constructed. For instance, the 
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FIG. 1 . Absolute errors (in degrees, right axis) for partial wave phase 
shifts with J < 4 (left axis) due to 7 different potentials fitting scat- 
tering data with ^ 2 /dof ~ 1 QJ-Q] as a function of the LAB energy 
(in MeV). 



triton binding energy obtained by Faddeev calculations is 
8.00,7.62,7.63,7.62,7.72 MeV for the CD Bonn d, Nijm- 
II, Reid93, Nijm-I and AVI 8 111 ill respectively. More recently, 
the covariant spectator model has produced the closest bind- 
ing energy 8.50 MeV to experiment precisely when the NN 
X 2 /dof is smallest. This yields in all B 3 = 7.85(34)MeV (exp. 
B 3 = 8.4820(l)MeV) i.e. AB 3 /3 = O.llMeV. Of course, in 
doing so even for the triton or the a-particle there is typically 
a flagrant need for three-body interactions which account for 
the missing lMeV and 4MeV to the binding energy respec- 
tively. On the other hand, the definition of the three- and 
higher-body interaction depends on the two body potential, 
so any uncertainty in the two-body interaction will carry over 
to the three-body interaction. Thus, even if we fix it say in 
the A = 3 system, there will always be a residual uncertainty 
in the A + 1 = 4 calculation. Thus, estimating the two-body 
uncertainty provides a lower bound on the total uncertainty if 
the correlations between the two- and three-body forces are 
ignored. The argument generalizes trivially to any A— body 
interactions and A + 1 —nuclei. 



Unfortunately, the procedure outlined above of using dif- 
ferent potentials stops beyond the A = 4 nucleus, due to com- 
putational and theoretical difficulties related to the form of 
the potential. From an ab initio viewpoint, only Monte Carlo 
calculations may go up to A = 10 when potentials are fixed 
to be r— dependent with a nonlocality in terms of the relative 
angular momentum operator 11121 1 1311 . For that good reason 
the Argonne potential saga has been constructed sticking to 
this representation and culminating in the AV18 potential [3], 
an updated version of AV14 containing charge-independence- 
breaking (CIB) terms and a complete electromagnetic inter- 
action and fitted directly to 4301 pp and np data from the 
1993 Nijmegen partial-wave analysis HI with the requested 
X 2 /dof ~ 1. The AVI 8 and AV18+UIX have become stan- 
dard Hamiltonians for ab initio calculations of light nuclei |Q] 
and dense matter IU4I1 . Note that even if statistical errors on 
potential parameters would have been estimated in Ref. J^, 
the question on the systematic errors remains. We motivate 
below an approximate method to address these issues. 

One of the outstanding features of the AV18-potential is the 
presence of the short distance core in the central part, Vc(r), 
for distances below a cole = 0.5fm, which demands sizeable 
and fine-tuned short distance correlations |6|. However, for a 
closed-shell nucleus one has schematically, 



(V2) A = 



MA- 1) r, 3 



d i rP 2 (r)V c (r), 



(1) 



where Piir) is the probability of finding two particles at a dis- 
tance r that turns out to be fairly independent of the particle 
number. In particular, for r < a core one has Vc(r) ^S> B/A and 
one is left with a two-body problem with A — 2 spectators; 
in the classically forbidden region an exponential suppres- 
sion, P2(r) ~ exp(— 2 f" core dr yj 2iiVc{r)) is expected semi- 
classically. Thus the contribution from the core is small, pre- 
cisely in the region where the NN force is not well determined 
from the PWA probing r > b m i n . For our error estimate we 
propose to side-step this core complication by introducing a 
coarse grained potential where the cancellation of the product 
P2{r)Vc{r) comes from a vanishing potential below a C ore- 

The previous argument was suggested long ago by Afnan 
and Tang [15] who realized that for A = 3,4 systems the 
relevant NN-scattering energies do not probe the core ex- 
pliticly. Soft core potentials, fitted to NN low partial waves 
up to £lab = 100MeV, provided reasonable binding ener- 
gies. Furthermore the hard core can be made into a soft 
core by introducing (linear-momentum) nonlocalities in terms 
of the kinetic energy operator by a unitary phase-preserving 
transformation [16]. Actually, this is the physics behind the 
so-called Vi ow k potentials based in the definition of an ef- 
fective truncated Hilbert space below a given cut-off A ~ 
^/Mnitijc 11171 Il8ll . Recently, we have shown how a similar 



idea can be implemented in coordinate space using a coarse 
grained potential lfl9l I20I1 . i.e. an average potential over a 
given wavelength resolution Ar ^ ^min 

; that means specify- 
ing the potential information in a finite number of points. The 
form of the potential is not important but calculations become 
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simple by taking delta-shells in the region below 3fm. For the 
partial wave 25+1 the potential reads 



compatible value, to ~ — Alt 1 / (M#A) J3- For nuclear matter 



Z i- L «=i 



r < r r 



(2) 



with /i the reduced pn-mass and r c — 3fm. In practice N < 5 
for any given partial wave. For r > 3fm we use the custom- 
ary CD OPE+electromagnetic interactions. The main novelty 
is a determination of the, so far, missing errors in the poten- 
tial parameters (in this case (A„)^,) from the uncertainties de- 
picted in Fig.Q]and corresponding to all % 2 /dof ~ 1 fits lUrjU. 
For instance, we found Il20ll by using the 3 5i, 3 L>i and E\ un- 
certainties of Fig. Q] that for the deuteron ABi H /2 = O.lMeV 
compared to the experimental AZ?2 H /2 = 0.005KeV. 

In our previous calculation [19], we showed how our ap- 
proach is competitive not only as a way of determining the 
phase shifts but also with more sophisticated approaches to 
Nuclear Structure lfl6ll . This was checked with oscillator wave 
functions in the case of 4 He, ie O and 40 Ca which reproduces 
experiment at the 20 — 30%-level provided the phase-shifts are 
fitted to about lOOMeV 11911 . This is a tolerable accuracy as 
we just intend to make a first estimate on the systematic un- 
certainties and then compute the change in the binding energy 
from the simple formulas, 

AB 3H = (AV 2 ) 3H = 3<1*|I (AV 1So +AV 3Si ) \ls) , (3) 

Afi 4H e = (AV 2 >4 H e = 6<1*| | (aVi S() + AV, Sl ) \1s) , (4) 

where is the Harmonic oscillator relative wave func- 
tion with the corresponding b— oscillator parameter repro- 
ducing the physical charge radius. The numbers in front are 
Talmi-Moshinsky coefficients and correspond in this particu- 
lar case to the number of pairs interacting through a relative 
s-wave. Errors are computed by adding individual contribu- 
tions {AXn) 1 , ,, from Eq. (O in quadrature. By propagating the 
PWA errors in Eq. ® we find AB 3 /3 = 0.07 - 0.085MeV 
depending on the fitting cut-off LAB energy, 100-350 MeV 
respectively, in good agreement with the Faddeev estimates 
given above. For the a— particle Eq. (0]l yields AB4/4 = 
0.10 - 0.13MeV. Along the lines of Ref. we also find 
ABi6 /16 = 0.26MeV and AS 4 o Ca /40 = 0.32MeV. 

A simple estimate on the impact of errors due to the two 
body interaction uncertainty can be done using Skyrme effec- 
tive interactions (for a review see [22]) 



AB 3 f , , N , 

~a = u At ° J d x P( x ¥> 



(5) 



where we get t = (w/m)E„»^(Ai Soj „ + h Su „). Using the 
two-body interaction of Eq. (O and propagating errors we 
get fo = — 0.92(l)GeVfm 3 . As a check on the size of fo 
we note that from a fit to the equation of state used by the 
Trento group 1E3I1 at low densities we obtain a value fo ~ 
-0.9(l)GeVfm 3 , whereas a coarse graining of NN interac- 
tions in CM momentum space down to A ~ 0.3GeV gives a 



at saturation, po = 0.17fm , our Afo = lOMeVfm implies 

(6) 



AB 3 

— = -AfoPo = 0.6MeV . 
A 8 



We may implement finite size effects by using a Fermi-type 
shape for the matter density p (r) =C/(1 +e^ R ^ a ) with R = 
roA3 and ro = l.lfm and a = 0.7fm and normalized to the 
total number of particles A = / d 3 xp (x) we get 



Afi A /A = 0.1-0.4MeV, 
which depends on the value of A for 4 < A < 208. 



(7) 
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FIG. 2. AVI 8 phase shifts in the 1 Sq (upper panel) and 3 5i (lower 
panel) partial waves as a function of the CM momentum when the pa- 
rameters are varied as specified in Ref. 12111 . We also depict the band 
corresponding to the spread of values obtained with 7 high quality 
potentials |I|43] containing One-Pion-Exchange (OPE) and Charge 
Dependence (CD) tail and fit scattering data with ^ 2 /dof ~ 1. 

One may reasonably doubt that the core effects and the cor- 
responding short distance correlations can be reliably moni- 
tored by a simple shell model calculation as given by Eq. (O 
and Eq. © for 3 H and 4 He respectively. We show now that 
for the purpose of error estimate this is however not so. In 
order to check this we take advantage of a recent analysis us- 
ing the AVI 8 -potential 112 in where the dependence of nuclear 
binding on the potential parameters is analyzed in much detail 
for mass numbers A = 2, ... 8. Relative 1% variations of two- 
body potential parameters are considered yielding changes in 
binding energies in the range of 0.1 — 20MeV for 4 He. This 
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Parameter 


e v 


e v (AVI 8) 


2 H 


2 H (AVI 8) 


3 H(100) 


3 H (350) 


3 H (AV18) 


4 He (100) 


4 He (350) 


4 He (AVI 8) 


m N 


-0.0445 


-0.0444 


-0.3997 


-0.3960 


-1.0030 


-1.1606 


-0.9797 


-1.6820 


-1.9973 


-2.1945 


mm + 8n 


-0.1602 


-0.1612 


-0.7986 


-0.7928 


-2.7817 


-3.8063 


-2.6560 


-5.2395 


-7.2887 


-6.5230 


8a 


0.1350 


0.1359 


0.4641 


0.4609 


2.0690 


3.0778 


1.9504 


4.1380 


6.1555 


5.0380 


m K (OPE) 


-0.0051 


-0.0051 


0.0289 


0.0623 


0.0140 


-0.1676 


0.0730 


0.0280 


-0.1461 


0.1045 


m n (+TPE-S) 


0.0706 


0.0705 


0.2665 


0.2981 


1.1205 


1.1770 


1.1055 


2.2409 


2.3540 


2.8105 


m„ (+TPE-L) 






0.2683 


0.2999 


1.1229 


1.1812 


1.1025 


2.2459 


2.3623 


2.7830 


my 


-0.4757 


-0.5079 


-1.8692 


-1.8571 


-7.9517 


-10.3458 


-7.6743 


-15.9034 


-20.6916 


-19.9870 



TABLE I. Changes of ground state energies of lightest nuclei A = 2,3,4 in MeV for the coarse-grained potentials with oscillator shell model 
when 1 Sq and 3 5i phase shifts are fitted up to LAB energies £Lab = 100MeV and £Lab = 350MeV used in this work compared to the ab 
initio Monte Carlo calculations |21] (labeled as AV18 here and using their notation) . The changes correspond to the variation of the total 
energy when the mass parameter is either increased or decreased by 1%. For instance the label m^ corresponds to B(l.Olmjv) — fi(0.99/n#). 



is a much larger range than our estimated errors. We note 
here that the largest sensitivity is on the short range potential, 
whereas the pion mass variation in the OPE piece yields a tiny 
effect. We may test our strategy by proceeding as follows. 
For any variation of the AV18-potential parameters there is a 
corresponding change in the phase shifts. In Fig. [2] we show 
as an illustration the changes in the most important 'So and 
3 5i waves due to several changes in parameters as explained 
in Ref. 12111 . As we see these changes are indeed larger than 
the systematic errors depicted in Fig.Q] that is enough for our 
purposes. Given this variation we then readjust our coarse 
grained potential, Eq. (|2), as to reproduce such a change and 
then use Eqs. ( 13141 ). We show our results in table U depending 
on the fitted maximal £lab- The disagreement of our results 
with those of AVI 8 in the case of changing the pion mass in 
the OPE potential is not very important since the net result is 
rather small anyhow and the pion mass is well known. We 
note, however, a larger sensitivity of the AVI 8 potential with 
respect to the short distance variations which is the relevant as- 
pect for our error analysis. In other words, our coarse grained 
calculation reproduces ab initio Afi's at the 20% accuracy. 

We summarize our points. Nuclear Binding energies are a 
crucial test for the Nuclear Many Body Problem. While first 
principles calculations are hampered by computational diffi- 
culties, we note that nuclear force uncertainties may have a 
useful impact on these calculations. The present theoretical 
estimates are in the range AB/A ~ 0.1 — 0.4MeV exceeding 
two or three orders of magnitude the avalaible precision of ab 
initio Monte Carlo calculations achieved up to date for A < 10. 
The prospective of increasing the particle number keeping the 
computational cost provides strong reasons to go beyond NN 
uncertainties as done here and to consider also errors in 3N 
and 4N forces. Finally, one should keep in mind that agree- 
ment between theory and experiment could also be declared 
even when the theory is less precise than the experiment, as it 
so frequently happens in Nuclear Physics, provided of course 
both error bands overlapp. 
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